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Abstract 

We present a diagrammatic analysis of baryons in the 1/A expansion, where A is the 
number of QCD colors. We use this method to show that there are an infinite number 
of degenerate baryon states in the large- A limit. We also show that forward matrix 
elements of quark bilinear operators satisfy the static quark- model relations in this limit, 
and enumerate the corrections to these relations to all orders in 1/A. These results hold 
for any number of light quark flavors, and the methods used can be extended to arbitrary 
operators. Our results imply that for two flavors, the quark-model relations for the axial 
currents and magnetic moments get corrections of order 1/A 2 . For three or more flavors, 
the results are more complicated, and corrections are generically of order 1/A. We write 
an explicit effective lagrangian which can be used to carry out chiral perturbation theory 
calculations in the 1/A expansion. Finally, we compare our results to what is expected 
from a chiral constituent quark model. 
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1. Introduction 



It is well-known that QCD has a non-trivial limit as N ^ 00, where N is the number 
of colors [1]. This limit can be defined by writing the QCD covariant derivative as 



where is the gluon field, and letting N — > 00 holding g fixed. QCD in the large- N limit 
is believed to share many features with QCD at N = 3, including confinement and chiral 
symmetry breaking. (It can be shown that confinement implies chiral symmetry breaking 
in the large- N limit of QCD [2].) Therefore, an expansion in 1/N around the large- N limit 
of QCD can be expected to give valuable insights into non-perturbative QCD phenomena. 

The properties of baryons in the large- N limit were first discussed in detail by Witten 
[3], who argued that baryons can be described by Hartree-Fock equations in this limit. This 
approach gave a consistent picture of baryons, and showed how various physical quantities 
scale with N. Most of the subsequent work on baryons in the large- N limit was done in 
the context of the Skyrme model [4] , which has been argued to be exact as N — > 00 [5] . 
The baryon spectrum in the Skyrme model consists of an infinite number of degenerate 
states. It can also be shown that the axial current matrix elements obey the static quark- 
model relations in the large- N limit of the Skyrme model [6]. This result was extended to 
all operators in ref. [7]. Another approach was initiated by Gervais and Sakita [8], who 
showed that consistency of chiral perturbation theory with the large- N limit implies the 
quark- model relations for the axial current matrix elements in the large- N limit. These 
methods were recently rediscovered and extended by Dashen and Manohar [9] and Jenkins 
[10], who considered 1/N corrections in this approach. 

In this paper, we will use diagrammatic and group-theoretic methods to analyze ba- 
ryons in the 1/N expansion. Using these methods, we derive the quark-model relations 
directly from QCD in a way which makes clear the role of quarks as the constituents of 
baryons. These methods allow us to enumerate the corrections to these relations to all 
orders in the 1/N expansion for the baryon masses and matrix elements of quark bilinear 
operators for any number of light quark flavors. (The methods used can be easily extended 
to general matrix elements.) For 2 flavors, our results imply that the quark- model rela- 
tions for the axial currents and magnetic moments get corrections of order 1/N 2 . (For the 
non-singlet axial currents, this result was obtained in ref. [9].) For 3 or more flavors, the 
results are more complicated, since the baryon representations depend on N in this case. 
We find that the corrections to the quark-model relations are generically of order 1/N. 

Combining these results with current algebra ideas, we then write an effective la- 
grangian which can be used to carry out a simultaneous chiral and 1/N expansion for 
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processes involving baryons and soft pseudo Nambu-Goldstone bosons. This lagrangian 
contains couplings that grow with N, so that it is not manifest that its predictions are 
physically sensible. We have verified the consistency of the predictions of this lagrangian 
in some non-trivial cases (including all cases considered in refs. [8] [9] [10]), but we have 
not proved it in general. However, we argue that this effective lagrangian is the most 
general one with the right field content and symmetry structure in the large- N limit, and 
we therefore believe that it must be consistent on general grounds. 

We conclude with some remarks comparing our results to what is expected from a 
non-relativistic constituent quark model. 

2. iV-Counting Rules for Baryons 

In this section, we derive our main results on the structure of baryon masses and 
matrix elements in the 1/N expansion. We begin by giving a rough idea of the nature 
of our argument. Consider a baryon matrix element of the form (B'\TO\ ■ ■ -O n \B). The 
operators (D\ , . . . , O n can change the spin and flavor of the quarks in the initial baryon 
consistent with invariance under the flavor symmetry SU(Np). If the velocities of the 
initial and final baryon states are equal, then the spin and flavor changes must also be 
invariant under SU(2) V , the group of spatial rotations in the common baryon rest frame. 
It turns out that invariance under the combined SU(2) V x SU(Np) symmetry is especially 
restrictive when N is large for states with total spin J <C N: the amplitude to change the 
spin or flavor of a quark is suppressed by 1/N for such states. This combinatoric result is 
responsible for the appearance of the static quark-model relations in the large- N limit for 
forward baryon matrix elements. 

To make this argument precise, our strategy is to write a diagrammatic expansion for 
the quantities we wish to study, and then find properties of this expansion that hold to 
all orders in g. This method is not strictly rigorous, since it is not known whether the 
quantities in question can actually be obtained by "summing all diagrams." However, the 
results derived here depend only on the iV-counting of the diagrammatic expansion, which 
we believe to be reliable. 

We consider the QCD states 

\&0)=B e Cl a 1 a 1 A 1 a a N a N A N \ U /i \ Z > 

where a\, . . . , a^v = 1, . . . , Np are flavor indices, «i, . . . , cun =T, I are spin indices, and 
Ai,...,An = l,...,iV are color indices. The operators in eq. (2) create a quark 
with definite flavor, spin, and color, in a perturbative 1-quark state \4>). We take \4>) to 
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have average momentum zero and momentum uncertainty of order Aqcd- Our results are 
insensitive to the precise nature of the state \<f>). The important feature of the state \B ) 
for our analysis is that it has the right quantum numbers to be a 1-baryon state at rest: it 
is a color singlet, has unit baryon number, and has angular momentum and flavor quantum 
numbers determined by the tensor B. 

The tensor B in eq. (2) is symmetric under combined interchange of spin and flavor 
indices, a±ai <-> a2«2, etc. B may be thought of as a quark-model wave function describing 
states with spin J = |, . . . ,|iV (for N odd). The Young tableaux of the SU(N F ) repre- 
sentation with spin J is shown in fig. (la). We will assume that the lowest-lying baryon 
states in QCD have these quantum numbers for arbitrarily large N, with the masses of 
the baryons increasing with J. This assumption is supported by the Skyrme model, and 
of course is true for N = 3. 

Under time evolution, the state \Bq) projects onto the state \B), where \B) is the 
minimal energy eigenstate of the full hamiltonian with a nonzero overlap with \Bq)J 

U(0, -T)\Bo) — \B)(B\B }e-^ T + ■■■, 
(B \U(T, 0) ^{B G \B){B\e- lE * T + • • • , 

as T^— zoo, where the ellipses denote terms which are exponentially suppressed. (The 
analytic continuation is the usual Wick rotation, which is justified by the assumption that 
the only singularities of the time-evolution operator are those required by unitarity and 
crossing.) Eq. (3) allows us to write matrix elements involving the state \B) as a sum of 
diagrams with external states of the form \B$). In our approach the state \Bq) plays a role 
similar to that of the perturbative vacuum, where a formula like eq. (3) is used to derive a 
diagrammatic expansion for time-ordered products in the interacting vacuum. 

If we choose the tensor B in eq. (2) to have definite spin-flavor quantum numbers, then 
by the assumptions above, there is a 1-baryon state with the same quantum numbers as 
| Bo). We would like to identify \B) with this 1-baryon state, but we must ensure that there 
is no state with lower energy with the same quantum numbers. Since \B) has unit baryon 
number, the candidates for this state contain a single baryon and any number (possibly 
zero) of mesons, glueballs, and exotics. For example, if \B) has spin J, we can consider 
states consisting of a spin- (J — 1) baryon and a meson. Since mesons have masses which 
are order 1 in the 1/N expansion, this state will have more energy than the 1-baryon state 

t For notational simplicity, we work in finite volume; the spectrum of states is then 
discrete, and we take all states to be normalized to unity. If \B) is a 1-baryon state with 
zero momentum, the overlap (B\Bq) goes to zero in the infinite volume limit as l/y/V. 
Because all of our results are independent of (B\Bq), this is harmless. 
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provided that Mj — Mj_i — »0 as iV^oo. Generalizing these considerations, we see that 
the method we are using allows us to study those states which become degenerate with 
the lowest-lying baryon states (assumed to be the J = \ multiplet) in the large-iV limit.''' 
We already know from previous work that there are an infinite number of such states in 
the large- N limit; this result will also emerge self-consistently from our analysis. 

2.1. Masses 



In this subsection, we derive our main results on the baryon masses. We will work 
in the limit where the current quark mass differences and electromagnetism are turned 
off, so that the theory has an exact SU(Np) flavor symmetry. In this limit, the baryon 
spectrum consists of multiplets with spin J = |, . . . ,|iV with the baryons degenerate in 
each multiplet. (We will show later how to treat SU(Np) breaking perturbatively using 
chiral perturbation theory.) 

Consider the quantity 

Z=(B Q \e-^ T \B ), (4) 
where H is the QCD hamiltonian. Using eq. (3), we have 

Y In Z^-iE B , (5) 

as T — zoo. To define a diagrammatic expansion for Z, we write the hamiltonian as 
H = Ho + 5H. We choose Hq to be a 1-body operator such that \Bq) is an eigenstate; 
consequently, 5H will in general contain 1-body "interactions," but these do not affect the 
iV-counting of the expansion. Once again, our results are independent of the precise form 
of H and 5H. We then write 

Z = e- iE ° T (Bo\Ui{\T,-\T)\B Q ), (6) 

where Hq\Bq) = Eq\Bq), and 



to) = Texp 



-% [ dtSH^t) 

Jt 



H T (t) = e iftot SHe- iAot , (7) 



t One might worry about states containing a baryon and soft pions in the chiral limit. 
However, we can add SU (Np) -preserving current quark mass term to make the pions mas- 
sive and remove the degeneracy of these states with the one-baryon states. The smoothness 
of the chiral limit then allows us to extrapolate our results to zero current quark masses. 
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Eq. (6) can be expanded as a sum of Feynman diagrams using a simple adaptation of 
standard techniques. Each term in the expansion of eq. (6) in powers of 8H consists of a 
time-ordered product of operators. Using Wick's theorem, this can be written as a sum 
of all possible normal orderings and contractions. The terms in which all operators are 
contracted are just the usual "vacuum" graphs which sum to give the vacuum energy. 
However, because the state \Bo) contains particles, the normal-ordered terms can also 
contribute. We represent these contributions diagrammatically by associating a dashed 
line for each external quark creation or annihilation operator, as shown in fig. (2b, c). In 
each such term, it is understood that all creation and annihilation operators are normal 
ordered. 

In this diagrammatic expansion, the mass of the state \B) is given by 

M B = E B - E vac = J2cr{B \6^\B ), (8) 

r 

where E vac is the vacuum energy and is an r-body operator of the form 

/n( r ) — y-oiar-OrOrAt ? l b 1 (3 1 A 1 ~t t,b r /3 r A r (q\ 

w ~ uX b 1 /3 1 ---b r /3 r LL a 1 a 1 A 1 a a a r a r A r a • V*) 

The external creation and annihilation operators again create a quark with definite spin, 
flavor, and color, in the 1-quark state \<j>) used to define the state \Bq). (We can take \<f>) 
to be one of a complete basis of 1-quark states, and write all contributions in terms of 
creation and annihilation operators in this basis. Operators which annihilate quarks in 
1-quark states other than \<f>) then vanish on states of the form |£> ).) To see that eq. (9) 
is the most general form of a color-singlet r-body operator, note that operators involving 
color e symbols can always be reduced to this form by using the identity 

e Al - AN 6 Bl ... BN =detA, where A at = 6*' t (s, t = 1, . . . , N). (10) 

Since QCD conserves angular momentum and flavor, the operators in eq. (8) must 
be singlets under SU(2) V x SU(Np), where SU(2) V denotes the group of spatial rotations in 
the baryon rest frame. The classification of these operators is a group-theoretical problem 
which is solved in appendix A. Since the results concern only the spin and flavor structure, 
it is convenient to express the results in terms of a spin-flavor Fock space: we define 

\B) = B aiai --- aNaN <a 1 ---*l NOtN \0). (11) 

Because the baryon wavefunction is symmetric under a\a.\ <-> 0,20.2 etc., the creation 
operators at* are bosonic. We use the "curved bra-ket" notation to distinguish these states 
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from QCD states. Given any QCD operator O appearing in eq. (8), we can construct an 
operator in the spin-flavor Fock space such that {Bq\0\Bq) oc (B\0\B), so we can write 

M B = E B - E vac = J2cAB\0^\B) : (12) 

r 

More details are given in the appendix. 

The result of the classification is that an arbitrary r-body operator appearing in eq. 
(12) can be written as a linear combination of 

{ir 2 vY- (is) 

Here {1} is the quark number operator and {J 2 } is the usual angular momentum Casimir 
operator: 

{1}| J) = N\ J), {J 2 }\J) = J(J + 1)\J), (14) 

where | J) denotes any state \B) with total spin J. We emphasize that this result holds for 
any number of flavors Np > 2. This result can be understood heuristically by noting that 
{J 2 } is a Casimir operator which suffices to label the possible representations of \B), since 
the flavor representations are determined by J (see fig. (la)). Of course, we can eliminate 
the operator {1} using the relation {1} = N which holds on states with baryon number 1. 

One way to determine the N dependence of the coefficients in eq. (8) is to use the fact 
that the baryon mass is order N [3] . Because we assume that the baryon masses increase 
with J, Mg must be greater than the mass of the J = \ baryon multiplet, which is of 
order N. On the other hand, Mb can be no greater than the mass of the 1-baryon state 
with the same quantum numbers, which is also of order N. Therefore, Mb ~ N for all 
states |£>), and each term in eq. (8) must be at most of order N. This implies 

Mb = Mo + c s {B\{J 2 Y\B\ c s < + • • • , (15) 

s 

where M ~ iV and the ellipses denote terms higher order in 1/N. Since there seems to 
be no reason to suppose otherwise, we will assume that the inequalities in eq. (15) are 
saturated. 

Eq. (15) is our main result for the baryon masses. We can immediately see that the 
masses of states with J ~ 1 become degenerate in the large- A?" limit. This allows us to 
conclude that the states with J ~ 1 are in fact 1-baryon states, as discussed above. The 
fact that the baryon spectrum consists of an infinite tower of degenerate states in the 
large- N limit is a well-known result of the Skyrme model, and is also a consequence of 
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the methods of ref. [8]. We also see that the leading mass differences between states with 
J ~ 1 are given by 

Mj -Mj,=^ [J( J + 1) - J' ( J' + 1)] . (16) 

This result was derived recently in ref. [10] by demanding consistency of chiral perturbation 
theory in the large- N limit. We have extended these results by classifying the corrections 
to all orders in 1/N. A qualitative picture of the baryon spectrum for large values of N is 
shown in fig. (lb). 

The previous discussion relied on the results of ref. [3] , and does not make use of the 
full power of the diagrammatic approach. We now consider the diagrammatic expansion 
in more detail and show that eq. (15) can be derived entirely within this approach. 

We begin by showing that InZ can be expressed as a sum of connected diagrams. 
The expansion of Z clearly includes disconnected diagrams such as the ones shown in fig. 
(2c). A diagram with n disconnected components is proportional to T n , where T is the 
time extent. This can be most easily seen in position space where there is one overall 
time integration for each disconnected component ranging from —\T to |T. However, 
since InZ = —iEsT, we know that the disconnected contributions must exponentiate in 
some way. For the vacuum graphs (i.e. those with no external creation or annihilation 
operators), this exponentiation is very simple: it is a standard result in quantum field 
theory that the disconnected diagrams have the right combinatoric factors to combine into 
an exponential. The graphs with external creation and annihilation operators have the 
same combinatoric factors, so we can write 

InZ = -iE vac T + ln(B Q \:e A :\B ), (17) 

where E vac is the sum of all connected vacuum graphs, and A is the sum of connected 
graphs with external creation and annihilation operators. If \Bq) is a state with definite 
total spin, then it is an eigenstate of any SU(2) V x SU(Np) singlet operator of the form 
of eq. (9), so that 

ln(B \:e A :\B ) = (B \ln:e A :\B ) 

= (Bo\:A:\B ) +l(B \:A 2 : - (:A:) 2 \B ) + ■■■ (18) 

We see that the disconnected diagrams with external creation and annihilation operators 
do not exponentiate directly. 

However, it is not hard to see that the right-hand side of eq. (18) is a sum of connected 
diagrams of a new type. Consider the 0(A n ) terms in the expansion of \nZ. By anticom- 
muting the quark operators, each 0(A n ) term can be written in a standard form as (:A:) n + 
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terms with anticommutators. The relevant anticommutator is Snix—y) = {ip~(x), ifj + (y)}, 
where ip is the quark field and we use ± to indicate the creation/annihilation part of the 
field. Sh is a Green's function, and is denoted diagrammatically by a dashed line. (Since it 
appears only in internal lines, there can be no confusion with the external creation and an- 
nihilation operators.) Physically, Sh is a correction to the propagator required to properly 
describe the propagation of "holes" (the absence of a quark in the initial state) . A typical 
diagram involving hole propagators is shown in fig. (2d). We now say that a diagram in 
this standard form is "connected" if it is connected by either Feynman propagators or 
hole propagators. A graph with n connected components is proportional to T n , so only 
the connected graphs can contribute to InZ; the disconnected diagrams must cancel. It 
is easy to see explicitly that this happens in the 0(A 2 ) term in eq. (18). This expansion 
in terms of connected diagrams with hole propagators is a generalization of the connected 
cluster expansion for Hartree-Fock perturbation theory which is well known in many-body 
physics [11]. 

The Accounting of diagrams is now very simple: we associate one power of l/\/~N 
to each gluon vertex and one power of N to each internal color loop. (Color loops are 
most easily counted using the double-line notation [1].) In this way, one can see that the 
diagrammatic expansion for the mass of the state \B) can be written 

M B = 5>(B|0 (r) |B), c r <j±-^ (19) 

r 

where is an r-body operator in the spin-flavor Fock space. Note that 

{1} = a\ a a aa = 1-body operator, (20) 
{J 2 } = {Jj}{Jj} = 2-body operator, (21) 

where {Jj} = a\ a { Jj)pa al3 . Eq. (19) therefore gives rise to the results for the masses 
quoted above. 

2.2. Matrix Elements 



We now extend the techniques described above to operator matrix elements. Using 
eq. (3), we can write^ 

mnrn n 1 n it M - «ggnMzMj°W im 



t We choose the phases of \B) and \B') so that (Bq\B) and (B' \B') are real. 
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Here, Oi_,...,O n are Heisenberg-picture operators, On, . . . ,Oj n are the corresponding 
interaction-picture operators, and we use the abbreviation Uj = t//(|T, — \T). The state 
has the same form as eq. (2), possibly with different spin and flavor quantum numbers. 
(In particular, it is defined in terms of the same 1-quark state \<f>) used to define |£>o).) 
The times ti, . . . , t n in eq. (22) must be lie within a fixed finite time interval as T — > oo in 
order for the states \B ) and \B' ) to project onto the corresponding interacting states \B) 
and \B'). 

We now consider the diagrammatic expansion for the right-hand side of eq. (22). 
The numerator can be written as the sum of diagrams with insertions of the operators 
0\, . . . , O n . The denominator can be written as a sum of diagrams with combinatoric 
factors obtained by expanding out the square roots. The numerator and denominator can 
be combined into a sum of diagrams involving hole propagators in a manner similar to 
that of eq. (18). The diagrammatic expansion therefore gives 

{B'\T6 1 {p 1 M)---6 n {p n ,t n )\B) = Y J Fz{^ (23) 

i 

where we have explicitly indicated the dependence on the 3-momenta of the operators. 
Each term in this series can be thought of as a form factor. The coefficient Fg contains all 
of the kinematic information of the form factor. 0[ is an 77-body operator of the form 
of eq. (9) (with no kinematic dependence) which gives the spin and flavor dependence of 
the form factor. 

We now argue that only connected diagrams contribute to the right-hand side of 
eq. (23). First, we note that we can discard any graph in which the insertions of the 
operators Oi, . . . O n are in different connected components by giving the operators definite 
3-momenta and restricting attention to "unexceptional" kinematic regions where no proper 
subset of the 3-momenta sum to zero. A graph with n connected components is therefore 
proportional to T n_1 , since the overall time integrations for the n—1 connected components 
which do not contain the operator insertions can range from —\T to |T; the connected 
component containing the operators has no overall time integral because the operators act 
at definite times. However, because (B'\TOi(ti) ■ ■ ■ O n (t n )\B) ~ T°, only the connected 
diagrams contribute. 

By counting powers of N in connected graphs, we find that the coefficient of an r^-body 
operator in eq. (23) is 

where L is the minimum number of quark loops in graphs that contribute to the matrix 
element. These results are valid as long as the momenta of the operators are held fixed 
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as N ^oo. In this case, the initial and final baryon velocity are the same in the large- N 
limit. If we attempt to describe processes in which the baryon velocity changes, some of 
the momenta of the operators must become large with N. In this case, there is additional 
N dependence coming from the overlap integrals between the initial and final states and 
the matrix elements are suppressed by e~ N [3]. 

As was the case for the baryon masses, we can now obtain interesting results from 
a group-theoretic classification of the operators that appear on the right-hand side of eq. 
(23). In terms of the spin-flavor Fock space, we have 

(B'lm^, *!)... O n (p n , t n )\B) = Fiipuh, • • • ,Pn, t n )(B\o\ ri) \B). (25) 

e 

Note that the coefficients Ft may have SU(2) V indices coming from non-trivial dependence 
on the angles of the momenta pi, . . . ,p n , so On does not necessarily have the same quantum 
numbers as the product of operators on the left-hand side. 

The matrix elements of most interest are forward matrix elements of a single quark 
bilinear of the form ipTWip, where ip is the quark field, W is a flavor matrix, and V 
is a Dirac spinor matrix. For such matrix elements, the operators On transform under 
SU(2) V x SU(N F ) the same way as the tensor WgTp. Here, is the static projection of 
r, i.e. 

(7iT*)? = (4)?, (7o)£ = <$, (26) 

etc. As shown in appendix A, the most general operator with these quantum numbers is 
a combination of 

{WT}, tv(W){T}, {W}{T}, {WJjHJMT}, (27) 

multiplied by powers of {J 2 } and {1}. Here, 

{WT} = ai a W b a T^, (28) 

etc. This classification holds for any number of flavors N F > 2. The operator {WT} 
gives rise to the quark model predictions for the matrix elements, which obey the well- 
known SU{2Np) "spin-flavor" relations. The methods used in the appendix for classifying 
operators can be easily extended to arbitrary operators. 

The dependence of the coefficients of the operators in eq. (27) can be read off from 
the rule in eq. (24). The operator tr(TT){r} has contributions from diagrams involving at 
least one quark loop, and so has coefficient 1/N. The other operators do not require quark 
loops, and so have coefficients 1/A r_1 for an r-body operator. 
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For N F > 2, the N- dependence of the matrix elements is somewhat complicated due 
to the fact that the size of the SU(2N F ) representation with given J 2 grows with N (see 
fig la). We will therefore begin by discussing the case Np = 2, where the baryon states 
are characterized by / = J = . . .. In this case, direct calculation shows that 

TORI—W, (29) 

for trW = 0. Therefore, the last operator in eq. (27) is redundant. t (This result does not 
hold for Np > 2.) 

We now consider the N dependence of matrix elements of the operators in states with 
J ~ 1. If both W and V are equal to the unit matrix (e.g. for the flavor singlet scalar 
density), then the matrix element is clearly of order N. If one of the matrices W and V is 
the unit matrix and the other is a traceless matrix (e.g. for the flavor non-singlet vector 
charge), the matrix element is of order 1, since {WT} is an SU(2) V x SU(Np) generator 
in this case. The only case which is not immediately obvious is when both W and V are 
traceless, corresponding to the flavor non-singlet axial current. Direct calculation shows 
that the matrix element is of order N in this case. 

These results, along with the order of the leading corrections to the large- N limit, are 
shown in table 1. Note that in every case, the quark- model result corresponding to the 
operator {WT} is exact in the large- N limit, and the corrections start at order 1/N 2 . 





w,r = i 


tiw = o, r = l 


W = 1, trT = 


tr W, tr T = 


{WT} 


N 


i 


1 


N 


tr(W){T}/N 


1* 





1/N* 





{W}{T}/N + h.c. 


N* 


i* 


V 


1/N 


{WT}{ J 2 }/N 2 + h.c. 


1/N 


1/N 2 


1/N 2 


1/N 



Table 1. Leading corrections to quark bilinear matrix elements for 2 light 
quark flavors. The corrections marked with a * are proportional to the 
lowest-order values. 

We now consider briefly the extension of these results to the case N F = 3. In this 
case, the SU(Nf) representations corresponding to spin J ~ 1 contain ~ iV states. We 

t We thank A. V. Manohar for pointing out this redundancy, which was missed in the 
initial version of this work. We also thank A. Cohen for discussions on this point, and on 
the A^-dependence of amplitudes involving baryons with finite strangeness in the 3-flavor 
case. 
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identify the physical baryon states with the large- N states that have the same isospin and 
strangeness. The iV-dependence of the matrix elements is now more complicated due to 
the dependence on the strangeness S of the baryon states. For example, the AS = 1 
semileptonic hyperon decays proceed through matrix elements of the QCD operator 

<?AS=1 = ^7 m *V>l, (30) 

where 
















|. x-| 














U 





0/ 



(31) 



In this case, the corresponding 1-body operator in the spin-flavor Fock space is {XJj}. 
This has matrix elements of order y/\S\N between states with strangeness S, since it 
removes a quark in a state with occupation number ~ S and places it into a state with 
occupation number ~ N. The computation of these matrix elements is rather complicated, 
and must be handled on a case-by-case basis. Generically, we find that corrections to the 
quark-model relations which depend on SU(3) f symmetry (and not only isospin symmetry 
alone) are of order 1/N. These issues will be discussed more fully in a future work, where 
we will also consider in detail whether the 1/N expansion is quantitatively relevant for 
N = 3. 

We close this section with some general observations on these results. The quark- 
model relations we have derived in the large- N limit are often said to be the consequence 
of an SU(2Np) spin-flavor symmetry. It is therefore natural to ask whether there is any 
sense in which QCD posseses such a symmetry in the large- N limit. Because any SU (2Np) 
symmetry should contain the SU(Np) flavor symmetry of the light quarks, all particles 
with SU(Np) quantum numbers must transform non-trivially under SU(2Np). t(This is 
in contrast with heavy-quark symmetry, which commutes with the chiral symmetry of the 
light quark flavors.) However, it is clear that the observed mesons (for example) do not lie 
in SU(6) multiplets. Because of this difficulty, we do not know of any sense in which the 
quark- model relations are the manifestation of an extra symmetry in the large- N limit. 

Finally, we note that it may seem paradoxical that "non-relativistic" relations such 
as the quark-model relations can arise in a relativistic system of quarks and gluons such 
as the one we have been studying. However, the quark-model relations we have derived 

t If all particles with SU(Np) quantum numbers decoupled from the baryons in the 
large-iV limit, then the SU(2N F ) symmetry could arise as a superselection rule in this 
limit. However, this does not happen. For example, the pion-nucleon scattering amplitude 
is order 1 in the 1/N expansion (see eq. (46) below). 
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are non-relativistic only in the sense that Lorentz invariance implies that they can hold 
in at most one reference frame. In our approach, this frame is picked out by the baryon's 
velocity, and the resulting picture is completely Lorentz covariant. 
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3. Chiral Perturbation Theory 



Using the results for operator matrix elements derived above, we can write down an 
effective lagrangian which describes the interactions of baryons with soft pions in the 1/N 
expansion. 

We begin by discussing the baryons. Because the baryon mass is of order N, we can 
describe the baryons using a heavy-particle effective field theory [12]. We write the baryon 
momentum as 

P = M v + k, (32) 

where M ~ N is a baryon mass and v is a 4- velocity (v 2 = 1) which defines the baryon rest 
frame. We then write an effective field theory in terms of baryon fields whose momentum 
modes are the residual momenta k. The reason that this effective theory is useful is that 
all kinematic dependence on M is removed. 

The Lorentz structure also simplifies considerably in a heavy-particle effective theory. 
All Lorentz invariants in the spin-| representation can be constructed from e^ uXp , and 
the spin matrix 

s m = 1( 7 m_^) 75 . (33) 

We find it convenient to simply work in the rest frame defined by v and use a non-relativistic 
notation. 

To describe the baryons with spin . . ., we use the baryon field B aiai "' aNaN (x). 
(In a relativistic notation, these fields would satisfy N constraints of the form P + B = 0, 
where the particle projection operator P + = |(1 +^) can act on any of the spinor indices.) 
The heavy field B is defined using a common mass M for all of the baryons, which we 
choose to be the mass of the J = \ multiplet in the SU(N F ) limit. With this choice, M 
is close to the physical mass of baryons with J ~ 1. 

The large- N relations derived in the previous section hold only for baryons with J ~ 1, 
while the effective lagrangian we are describing also includes baryon fields with J ~ N. 
Because these states have masses M + O(N), their effects on the dynamics of the states 
with J ~ 1 will be local (on the scale of external momenta which are order 1 in the 1/N 
expansion) and suppressed by powers of 1/N. (We do not want to consider an effective 
theory with these states integrated out, since integrating out a subset of the baryon states 
violates the spin-flavor relations explicitly.) Furthermore, baryons with J ~ N will only 
appear in intermediate states if we consider the insertion of operators with spin ~ iV or if 
we carry consider diagrams involving ~ iV vertices. However, it is still important to know 
whether these contributions, when they appear, signal a breakdown of the 1/N expansion 
as defined by this effective lagrangian. If the contributions from baryons with J ~ N have 
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the same form as the tree-level terms, then these contributions are harmless. While we 
believe this to be the case, we have no formal proof. We will see below that the consistency 
of the effective lagrangian we are describing is non-trivial even at low orders in the 1/N 
expansion. 

We now discuss the light mesons. The main new feature of chiral dynamics in the 
large- N limit is that the effects of the axial anomaly are subleading in 1/N, so that the 
pattern of spontaneous symmetry breakdown in the large- A?" limit is [13] [2] 

SU(N F ) L x U(1) L x SU(N F ) R x U(1) R — > SU(N F ) L+R x U(1) L+R . (34) 

Therefore, in addition to the usual Np — 1 Nambu-Goldstone bosons ( "pions" ) , there is an 
extra Nambu-Goldstone boson in the large- N limit, which can be identified with the rj'. 
At higher orders in 1/N, £7(1) l-r is broken explicitly. In this section, we restrict ourselves 
to making some remarks on chiral perturbation theory which do not involve the rj' in an 
essential way, so we will ignore the presence of the light rj' in the following discussion. 

The Np — 1 pions are described by a field 

e(x) = e m ^/^, (35) 

which is taken to transform under SU(N F ) L x SU(N F ) R as 

£ L£L/t = UiR\ (36) 

where this equation implicitly defines U as a function of L, R, and £. The effective 
lagrangian is most conveniently written in terms of the hermitian fields 

V, = l - (£«9 M £t + £ta M £) , A[x = % _ (^t _ e t^) j (37) 

transforming as 

V^UV^ -id^Utf, A^UA^UK (38) 

With these definitions, it is easy to write down the effective lagrangian. In order to 
keep track of the tensor algebra, it is convenient to once again utilize the spin-flavor Fock 
space. For example, we write the heavy field B in terms of 

\B) = B aiai --- aNaN < ai ---al NaN \0). (39) 

The chiral covariant derivative acting on baryon fields can then be written as 

V fl \B) = (d»-i{V fl })\B), (40) 
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and the first few terms in the effective lagrangian involving the baryon fields are 

£eff = (B\iV \B) + g(B\{A ■ s}\B) + c(B\{m}\B) + • • • , (41) 
where m is the quark mass matrix spurion, and 

{A ■ s} = {s%{Aj)>La b ^ {m} = mla\ a a b «. (42) 

To find the N dependence of the coefficients g and c in eq. (41), we use the fact that 
the terms in the effective lagrangian can be related to matrix elements of QCD operators. 
The term proportional to g gives rise to a contribution to the matrix element of the axial 
current, while the term term proportional to c gives rise to a contribution to the matrix 
element of the scalar density. Comparing with eq. (24), we see that g and c are both of 
order 1 in the 1/N expansion. In general, it is not hard to see that a general term in the 
effective lagrangian involving the baryon fields will have the form 

— ^ (B | r-body operator|£), (43) 

where the r-body operator in the spin-flavor Fock space is constructed out of the fields 
and spurions in the effective lagrangian. This form ensures that the iV-counting of the 
effective lagrangian is the same as that found in the previous section. 

We now consider the question of the consistency of this effective lagrangian. Suppose 
we form a graph from two terms in the effective lagrangian, corresponding to an r±- and 
an r2-body operator. The graph will therefore have the form l/N ri+T2 ~ 2 times an (n + 
r2)-body operator, apparently violating the iV-counting of the previous section. Such 
contributions can only be consistent if the (ri + r 2 )-body operator has a special form, so 
that it can be written as a fewer-body operator. We will show that this is in fact what 
happens in some examples, but we have not been able to prove the consistency of this 
effective lagrangian in general. 

We first consider baryon-pion scattering at low energies. For simplicity, we work in 
the chiral limit. The three graphs that contribute are shown in fig. (3). The graph in fig. 
(3c) gives a contribution of order 1/N to the amplitude, ^ while the graphs in fig. (3a,b) 
lead to an amplitude [8] [9] 

A(Btt? -> £Vf ) oc 4^ (B'\[X Ai , X Bk ]\B) + ■■■, (44) 

t We work with non-relativistically normalized states, so that there is no dependence 
on the mass in the normalization of amplitudes. 
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where 

X* A = {\ A ) a h {J%a\ a a b ^ (45) 

is the axial current in the spin-flavor Fock space corresponding to the pion -k a . Because 
the initial and final baryons are on shell, their residual momenta are 0(1/N), and can be 
neglected. Since {B'\X A ^\B) ~ A" and f % ~ VN [3], one naively expects A ~ N, which 
would violate unitarity for any fixed kinematics as N^oo. However the commutator of 
two axial currents is a vector current, so for baryons with J ~ 1 the term shown in eq. 
(44) is of order 1/N. If we include the mass splitting between baryons with different spins, 
we obtain a contribution to the amplitude proportional to 

^ (B'\X Aj {J 2 }X Bk + X Bk {J 2 }X A] \B) ~ 1, (46) 

J IT 

and so the full amplitude is actually of order 1 in the 1/N expansion. t 

Similarly, we can consider scattering processes involving more pions. The graphs 
shown in fig. (4) give rise to a scattering amplitude 



A(Bn A ^B>nM) oc [f 2 (B'\[X A \ [X B \X C %B) 

-q1(B , \[X Bk ,[X Ci ,X A *]]\B)) + 



(47) 



The naive estimate of this term is A ~ A^ 3 / 2 , but the double commutators lead to matrix 
elements of the form (B'\X\B) ~ N, so that in fact A ~ l/v^/V [9]. The finiteness of 
this amplitude in the large- N limit was used in ref. [9] to argue that there are no 1/N 
corrections to the axial current matrix element. It is easy to check that the corrections to 
the axial current that we have enumerated give rise to corrections to this amplitude which 
are consistent in the large- N limit. 

In both of the cases considered above, products of operators that naively violate the 
iV-counting rules had a commutator structure that allowed them to be expressed in terms 
of operators with fewer creation and annihilation operators. As argued above, a mechanism 
of this sort is required for consistency in the large-A^ limit. This commutator structure 
has its origin in the antisymmetry of tree amplitudes under crossing. For loop amplitudes, 
the "ze" part of the baryon propagator ruins this antisymmetry, and amplitudes are no 
longer proportional to commutators. For example, if the pions have a common mass m^, 
the baryon masses have a 1-loop contribution [10] 

AM * K TTTk im Aj X Aj \B) oc \N 2 + J (J + 1)] (48) 

loir f 4 l^ti 



t We thank A. V. Manohar for pointing this out to us. 
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for a baryon with spin J (fig. (5)). The leading term grows with N, but it is a harmless 
constant mass shift consistent with the iV-counting results given in table 1. The mass 
difference between baryons with different spin is of the form derived in section 1. 

We have verified that all of the examples considered in refs. [8] [9] [10] give consistent 
results using our approach. However, other than the remarks made above, we do not have 
any general insight into the cancellation mechanism at work. However, we believe that 
the effective lagrangian we have written down is the most general one compatible with the 
symmetries and iV-counting structure of the large- A?" limit. We believe that any effective 
lagrangian which is defined by its field content and symmetry structure is consistent as 
long as all allowed terms are included. (This "theorem" is advocated by Weinberg as 
the foundation of the method of effective lagrangians [14].) If this is correct, the effective 
lagrangian we have described provides an explicit and consistent description of pion-baryon 
interactions at low energies in the 1/N expansion. 

4. The Non-relativistic Constituent Quark Model 

Traditionally, SU(6) relations have been justified by appealing to a non-relativistic 
constituent quark model (NRCQM) in which the SU(6) symmetry arises because the 
constituent quarks are heavy. In this section, we make some brief remarks comparing our 
results to what is expected from such a model. We include this discussion to point out 
that the pattern of corrections to SU(6) relations we have derived above is very different 
from what is expected from a NRCQM. 

The iV-counting of the NRCQM is identical to that of QCD, and so 5 , L r (6)-violating 
amplitudes will be suppressed by powers of 1/N in the NRCQM. However, when we discuss 
the NRCQM below, we will not treat 1/N as a small parameter, since the point of our 
discussion is whether the SU(6) symmetry can be explained solely in terms of the heaviness 
of the constituent quarks. 

The NRCQM can be formulated by assuming that below the QCD chiral symmetry 
breaking scale but above the confinement scale, QCD dynamics can be approximated by 
an effective lagrangian containing constituent quarks and weakly-interacting gluons [15]. 
The constituent quarks are to be viewed as composite objects whose properties are the 
result of non-perturbative QCD dynamics at the chiral symmetry breaking scale. The 
most important difference between the constituent quarks and the "current" quarks that 
appear in the fundamental QCD lagrangian is that the constituent quarks have a flavor- 
independent mass Mq ~ Aqcd due to chiral symmetry breaking. Because chiral symmetry 
is broken spontaneously, the theory also contains 8 "pions" in the chiral limit. 

In this model, it is easy to see how the SU(6) spin-flavor symmetry arises. If we 
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assume that the effective gluon coupling a s is weak, then the constituent quarks inside a 
baryon will be nonrelativistic with binding energy E B ~ o? s Mq <C Mq. The interactions 
of the constituent quarks at energies of order Eb therefore exhibit an SU(6) spin-flavor 
symmetry exactly analogous to the spin-flavor symmetry discussed recently for b and c 
quarks [16]. This SU(6) symmetry will give rise to the same relations for baryon masses 
and matrix elements derived above in the large- N limit of QCD. 

However, it should be emphasized that the picture of QCD dynamics contained in the 
NRCQM has some serious drawbacks when taken seriously. First of all, if the constituent 
quarks are heavy and weakly interacting, there should be weakly bound states with the 
same quantum numbers as the pions with masses ~ | of a baryon mass. Furthermore, 
because of the SU(6) symmetry, these states should lie in an SU(6) multiplet with a nonet 
of vector mesons. Identifying these states with observed mesons is rather problematic. 
Another difficulty with the constituent quark model is that it leads us to expect hadron 
form factors to deviate significantly from point-like behavior for momentum transfers of 
order Eb <C Mq ~ 300 MeV. Instead, we find that baryon form factors vary on a scale 
of order m p ~ 770 MeV. Finally, the gluons which couple the constituent quarks confine 
at some scale below Eb (since the gluons are supposed to be weakly interacting at that 
scale), so we expect glueball states with mass <C Mq. It is not clear whether any of these 
objections is fatal, but we want to emphasize that the validity of the NRCQM should not 
be taken for granted, despite its successes. 

Leaving these problems aside, we want to compare the expected pattern of corrections 
to the SU(6) relations from the NRCQM to what we found for the 1/N corrections in QCD. 
This is easy to do if we write an effective lagrangian for the constituent quarks treated as 
heavy fields: 

£nrcqm = Qiv ■ DQ + g Q QA ■ sQ + -^-Qv^s v G^Q + 0(1/Mj). (49) 

ZMq 

Here Q is the constituent quark field, is the gluon covariant derivative, and G^ v = 
^iiv\ P G Xp is the dual of the gluon field strength. The SU(6) symmetry is violated by the 
pion couplings and by the color magnetic moment term. Pion loops give rise to SU(6)- 
breaking corrections to matrix elements which vanish at zero momentum in the chiral 
limit. When the current quark masses are turned on, the pion loops give 5 , L r (6)-breaking 
corrections at zero momentum proportional to powers of the current quark masses. The 
I/Mq terms give rise to 5 , L r (6)-breaking corrections at zero momentum. These corrections 
will be proportional to powers of the SU (6) -breaking spurion s^/Mq; only even powers 
of this spurion can appear in corrections by parity invariance. Thus, all of the SU (6) 
relations for matrix elements and baryon masses are expected to have 1/Mq corrections 
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in the NRCQM. This pattern of corrections is very different from that found in the 1/N 
expansion of QCD, where quark-model relations can have corrections which are either 1/N 
or 1/N 2 . 

5. Conclusions 

In this paper we have developed a diagrammatic expansion for baryons in the large- N 
limit. This method gives a simple recipe for determining the iV-dependence of physical 
quantities. Using these techniques, we have shown that forward matrix elements of quark 
bilinears obey SU(2Np) spin-flavor relations. Our approach makes a direct connection 
between the emergence of these relations and the spin-flavor properties of the baryon 
wavefunctions, which we feel is very attractive. 

On a more practical level, our methods allow the enumeration of the 1/N corrections 
to the large- N results to all orders in 1/N, completing previous partial results on the clas- 
sification of 1/N corrections. We have also written down explicitly an effective lagrangian 
which describes the interactions of baryons and pions at low energies in the 1/N expansion. 
The consistency of this lagrangian in the large- N limit is not manifest, because some of 
the couplings grow with N. While we have not been able to prove consistency in general, 
we have verified it in a number of cases and given general arguments indicating that this 
effective lagrangian must in fact be correct. 

After completing this paper, we received ref. [17], which considers baryon masses and 
matrix elements using the Hartree-Fock approach of ref. [3]. We also received ref. [18] 
which develops further the approach of refs. [8] [9]. 
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Appendix A. Classification of Operators 

In this appendix, we consider the classification of few-body operators with given trans- 
formation properties under SU(2) V x SU(Np). Since only the spin and flavor quantum 
numbers are relevant for our discussion, it is useful to use a notation in which only spin 
and flavor indices appear. To do this, we note that the fact that the states \Bq) are an- 
tisymmetrized in color allows us to treat the quarks as colorless bosons. To make this 
precise, we make use of the identity 

JV 

«U^l^o) = ^ lQ ™^e^-^^<^ al iaiAi ... a\ lAs ... a[ NaNAN \0). (50) 

s=l 1 s N 

Exactly the same algebraic relation is satisfied by bosonic creation and annihilation oper- 
ators carrying no color. Specifically, we define 

\B) = B aiai --- aNaN al ai ---al NOtN \0), (51) 

where the ct^'s are bosonic creation operators which act on a spin-flavor Fock space. We 
use the "curved bra-ket" notation to distinguish these states from QCD states. We then 
have 

al 1 a^\B)=B^--- a ^Y j 5 b a3 8l < Qi • • • ■ • • a^jO). (52) 

s=l 1 s N 

Comparison of eqs. (50) and (52) shows that 

(B \6^\B )cx{B\O^\B), (53) 
where we use the caret notation to denote QCD operators, and 

Our job is now to classify the general r-body operators with specified SU(2) v xSU(Nf) 
transformation properties. First note that the most general r-body operator can be written 

as 

o {r) = Kt'ZT^ ■ ■ ■ «l Q y 1/31 • • • ^ + ■■■, (55) 

where the omitted terms arise from the re-ordering of the creation and annihilation op- 
erators. These omitted terms are irrelevant for the classification of the operators, since 
they are s-body operators, with s < r. We can therefore classify the operators inductively, 
starting with 1-body operators. We will make frequent use of this simplification in this 
appendix. 
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With the ordering chosen above, the creation operators project out the part of X which 
is symmetric under simultaneous interchange of spin and flavor indices a\OL\ <-> a^cti, etc. 
For notational simplicity, we will consider the case of 3 flavors, although the results are 
valid for an arbitrary number of flavors. 

6.1. Singlets 



In this subsection, we classify the SU(2) V x 577(3) singlet r-body operators. We 
proceed inductively in r. It is clear that the only 1-body singlet operator is a\ OL a aa . 

Now suppose that we have classified all of the singlet operators up to the (r — l)-body 
operators, and consider the r-body operators. The tensor X in eq. (55) must be a linear 
combination of products of the 577(3) and SU(2) V invariant tensors 8%, e a6c , 5p, and e a ^. 
First we note that the spin and flavor e symbols always appear in pairs and can thus be 
replaced by Kronecker <5's using eq. (10). Therefore, the tensor X in eq. (55) can be written 
as a linear combination of terms of the form 

x hT'Tt = C • • • KV'b^ (56) 

0\p\---O r p r Oi O r pi"'Pr v ' 

where we have used the symmetry properties of the creation and annihilation operators 
to put the flavor indices in canonical order. The tensor Y is an SU(2) V singlet, and can 
therefore be written as a linear combination of tensors of the form 

where a is a permutation of 1, . . . , r. (We have used up our freedom to re-order the indices 
in writing eq. (56).) 

We can rewrite SU(2) V tensors of the form eq. (57) by writing the permutation a as a 
product of simple interchanges. In this way, Y can be written as a product of the tensors 

We can then use the identity 

(WMj)t 2 =^tj:-¥p1 6 t,} ( 59 ) 

to write Y as a linear combination of products of SU(2) V generators and the identity 
matrix. Finally, some simple SU{2) algebra gives us an expression for as a linear 
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combination of products of {Jj}'s and {l}'s (the number operator), which can be reduced 
to the form {l} r_2s { J 2 } s , where we have used the notation 

{Jj} = aUj^a a ^ (60) 

and {J 2 } = {JjHJj}. 

To understand these last steps in detail, consider a simple example which can be easily 
extended to the general case: 

P1P2P3 P2 P3 Pi P1P2 72P3 V ' 

Using eq. (59), we can write the corresponding 3-body operator as 

= ^{JjHJkJjHJk} - 2{i}{ J 2 } + H 1 ) 3 + • • • • ( 62 ) 

The ellipses in eq. (62) indicate 1- and 2-body operators obtained by re-ordering the 
creation and annihilation operators. Using the identity 

JjJk — \&jk + (63) 

we obtain 

0(3) = -{!}{ J 2 } + 2ie jke {J j }{J k }{J e } + 1{1} 3 + • • • . (64) 
Finally, we use the identity 

ejkt{Jk}{Jt} = \tjki [{Jk}, {Jt}] = i{Jj} (65) 

to write 

0( 3 ) = -3{l}{J 2 }+i{l} 3 + ---. (66) 

These steps can be repeated for an arbitrary tensor Y of the form of eq. (57). After 
writing Y as a product of the S tensors defined in eq. (58), we use the identity eq. (59) to 
write as a linear combination of products of {1} and operators of the form { J • ■ ■ J}, 
where the indices on the J's are contracted using <%'s. Using eq. (63), we reduce this 
to products of {l}'s and {Jj}'s, contracted with Sj^s and e^'s. We then use eq. (10) 
to eliminate the ejfc/s in pairs in favor of Sj^s, and then apply eq. (65) to eliminate the 
last €jkt (if any). In this way, we obtain products of {l}'s and {Jjj's contracted with 
SjkS, which can be written as {l} r 2s ({J 2 }) s up to operators with fewer creation and 
annihilation operators. Thus, we have shown that the most general r-body operator is a 
linear combination of operators of the form 

{ir 2 v 2 r- (67) 
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6.2. Quark Bilinear s 



In this subsection, we classify the possible operators with the quantum numbers of 
the tensor W^T^, corresponding to the matrix elements of a quark bilinear. An arbitrary 
quark bilinear operator can be written as a linear combination of operators of this form. 
Our argument is very similar to the one presented above for singlets, and so we will be 
somewhat terse. 

We again proceed inductively. We begin by considering the case tr W = tr V = 0. In 
this case, the only effective 2-body operator with the right quantum numbers is 

= a\ a W^a h K (68) 

Now assume we have classified all operators up to (r — l)-body operators, and consider the 
r-body operators. A general r-body tensor of the form of eq. (55) with the right quantum 
numbers can be written as a linear combination of tensors of the form 

■y-a\a\---a r a r Tl/b' W^' V a ' a ' a i a i ' ' •o, r a. T ( PS\\ 

^b 1 f3 1 ---b r f3 r — VV a' L a' I b'f3'b 1 p 1 ---b r p r > l Dy J 

where Y is an SU(2) V x £77(3) invariant tensor. It must therefore consist of combinations 
of 8% : s, 5p% and e symbols for the groups SU(2) V and 577(3). However, since Y has an 
equal number of upper and lower indices for both groups, the e tensors must appear in 
pairs, and we can use eq. (10) to eliminate them in favor of combinations of 5's. The tensor 
in eq. (69) thus reduces to a linear combination of the form 

bi/3i—b r /3 r bi b 2 b r fit p 1 ...p t .../3 r ' ^ ' 

where s,t = l,---,r, and the caret denotes omission. Y is an SU(2) V invariant tensor, 
formed from a linear combination of products of #g's, where the lower indices are permuted 
relative to their canonical order. Going through the same steps as for the singlet operators, 
we can write this tensor in terms of S^s and Jj's with indices in canonical order. In this 
way, an arbitrary r-body operator can be written as a linear combination of operators of 
the form 

£>M = {WJ ■ ■ ■ JTJ ■ ■ ■ J}{J ...J}.-. {J... J}{1} ... {1} + ... , (71) 

where the J's are contracted using Sj^s, {WT} = a\ a W^T^a bl3 , and we have discarded 
all lower-body operators. 
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We now use the identity eq. (63) to reduce this to a product of {l}'s and {JjYs 
contracted with Sj^s and e^'s. If T = Jj, the operator is a linear combination of 



0«~{W,}{J}...{J}, (72) 

{W}{J j }{J}---{J}, (73) 

{WJ}{J j }{J}---{J}, (74) 

6j...{W}{J}---{J}, (75) 

e j ...{WJ}{J}---{J}. (76) 



The unlabeled J's in eqs. (72)-(76) are contracted using and e^. We have dropped 
all factors of the operator {1}, since {1} = N on the states we are interested in, and these 
factors therefore do not affect the iV-counting (see eq. (24)). 

We now eliminate the e symbols pairwise using the identity eq. (10). This leaves at 
most one e symbol. The identity eq. (65) can be used to eliminate the final e symbol (if 
any) from eqs. (72)-(75). The only non-trivial way an e can appear in eq. (76) is as a 
factor of 

e jk e{WJk}{Je} = i [{WJj}, {J 2 }] +■■■, (77) 

where the ellipses denote 1-body operators. However, the right-hand side of eq. (77) can 
be seen to have the opposite transformation property as {WT} under time reversal, so this 
operator cannot appear. 

In this way, we find that the most general r-body operator with the quantum numbers 
of WY is a linear combination of the operators 

{WT}, {W}{T}, {WJMJjHT}, (78) 

multiplied by an SU(2) V x SU(3) singlet operator. (We may need to add or subtract the 
hermitian conjugate from the operators in eq. (78) so that the effective operator has the 
same hermiticity properties as the corresponding QCD operator.) 

It remains only to consider the case where W = 1 and/or T = 1. Repeating the above 
argument with trivial modifications shows that eq. (78) gives the most general corrections 
in this case as well. 
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Figure Captions 



Fig. 1. (a) The Young tableaux describing the SU(2) V x SU(Np) representation of the 
baryon states with spin J. (b) The spectrum of baryon states in the large- A?" limit. The 
mass splitting AM between adjacent states with spin J ~ 1 is of order 1/iV, while for 
J ~ VN, AM ~ 1; our methods do not allow us to study baryon states with J ~ N. 

Fig. 2. Typical diagrams contributing to the diagrammatic expansion of Z in eq. (4). (a) 
A "vacuum" graph arising from the fully contracted part of the Wick expansion, (b) A 3- 
body operator arising from the normal-ordered terms in the Wick expansion. The external 
creation and annihilation operators are denoted by dashed lines. The solid internal lines 
denote the usual quark propagators, (c) A disconnected diagram with external annihilation 
and creation operators. The entire contribution is normal-ordered, (d) A typical diagram 
involving the "hole" propagator Sh(% — y) (denoted by an internal dashed line). 

Fig. 3. The leading contributions to the amplitude for the process Bti — > B'n. The 
diagrams (a) and (b) are proportional to axial current matrix elements (see eq. (44)). 
Diagram (c) is proportional to the matrix element of a vector current, and is of order 
for baryons with J ~ 1. 

Fig. 4. Contributions to the amplitude for Btv — > B'tttt which are proportional to matrix 
elements of three axial currents. 

Fig. 5. A 1-loop diagram contributing to the baryon masses. 
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